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We investigate cosmological density perturbations in a covariant and gauge-invariant formalism, 
incorporating relativistic causal thermodynamics to give a self-consistent description. The gradient 
of density inhomogeneities splits covariantly into a scalar part, equivalent to the usual density 
perturbations, a rotational vector part that is determined by the vorticity, and a tensor part that 
describes the shape. We give the evolution equations for these parts in the general dissipative case. 
- - - Causal thermodynamics gives evolution equations for viscous stress and heat flux, which are coupled 

, to the density perturbation equation and to the entropy and temperature perturbation equations. 

' We give the full coupled system in the general dissipative case, and simplify the system in certain 

I I cases. A companion paper uses the general formalism to analyze damping of density perturbations 

■ before last scattering. 
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! I. INTRODUCTION 

(N . 
> ■ 

QQ . The analysis of density perturbations in cosmological fluids is well established, particularly using Bardeen's gauge 
invariant formalism This formalism is inherently linear (i.e., it starts from the background and perturbs away 

from it) and non-local. An alternative approach, developed by Ellis and Bruni is covariant (and therefore local) 
C and readily incorporates nonlinear effects (since it starts from the real spacetime, not the background). We will use 

, this covariant and gauge-invariant formalism, in which the variables have a clear physical and geometric interpretation. 
Furthermore, the covariant approach is directly compatible with causal relativistic thermodynamics, as developed by 
Israel and Stewart 

O ■ Although dissipative terms representing viscosity and heat conduction have been formally incorporated into the 
^ equations in both approaches , most applications of the theory are restricted to the non-dissipative case - and even 
5-H . in this case, relativistic thermodynamics is usually not applied to analyze the behavior of the fluid self-consistently. 

This is is not a problem when studying the evolution of large-scale perturbations, which are unaffected by local physics 
^ , - although the generation of these perturbations, their initial evolution before leaving the Hubble radius, and their 
! ' final evolution after re-entering the Hubble radius, are governed by local physics. For small-scale perturbations, within 
rN \ the Hubble radius, a self-consistent analysis requires the application of thermodynamics]^ 

■ Here and in a companion paper we develop and apply a covariant and gauge- invariant analysis of density 
■ ■ ' perturbations that self-consistently incorporates relativistic causal thermodynamics. The general evolution equations 

governing density inhomogeneities are considered in Sec. II. Inhomogeneities are covariantly characterized by a scalar 
part, which represents the usual density perturbations, a vector part, which we show is determined by the vorticity, 
and a tensor part, which determines the shape of gravitational clustering. New evolution equations are derived for the 
vector and tensor parts, as well as for perturbations in the number density, entropy and temperature. We use the Gibbs 
equation to incorporate the temperature and entropy self-consistently, and we covariantly characterize different types 
of perturbation. In Sec. Ill, the viscous stress and heat flux that appear in the perturbation evolution equations are 
subject to thermodynamical transport equations, which then form a coupled system with the perturbation evolution 
equations. We define appropriate dissipative scalars to obtain a closed system of dynamical equations. The equations 



*email: maartens@sms.port.ac.uk 
^email: pepOulises. uab.es 

'^We are considering here the case of hydrodynamics. Dissipative effects on the microwave background have been self- 
consistently analyzed via numerical integration of the Boltzmann equation (see, e.g., Q). A covariant and gauge- invariant 
approach to the Boltzmann equation is developed in [M. 
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are simplified in the particular cases of entropy perturbations (non-dissipative), and when only one form of dissipation 
is present. 

The Israel-Stewart transport equations are under reasonable conditions causal and stable and thus provide a 
consistent relativistic description of local physical effects on small-scale perturbations. The thermodynamics of Eckart 
(and a similar alternative due to Landau and Lifshitz) is more established in the literature. However, in this theory 
the transport equations reduce from evolution equations to algebraic constraints on viscosity and heat flux, and as a 
result, the theory is non-causal (dissipative effects propagate at super- luminal speeds), and all its equilibrium states 
are unstable It can be argued that these pathologies only occur outside the hydrodynamic regime. But firstly, 
the stability problem persists in all situations, and secondly, it seems preferable to employ a theory with built-in 
causality and stability. Furthermore, the causal theory can deal with transient and short-wavelength effects, which 
are important in many cosmological and astrophysical situations (see, e.g., |[lT| - p^ ]). 

Applications to dissipative situations are treated in a companion paper where we analyze viscous damping of 
density perturbations before last scattering. This generalizes the results of Weinberg who used non-causal Eckart 
thermodynamics. 



II. COVARIANT APPROACH TO DISSIPATIVE PERTURBATIONS 

The covariant and gauge-invariant analysis of density perturbations is fully discussed in |^,^ . Here we present only 
the main points that are necessary for our purposes, before going further by deriving new evolution equations and 
incorporating causal thermodynamics. Our notation and conventions follow |^,|l6|,0, with some changes (see [p^). 

Given a covariantly defined fluid four- velocity u"^ (see the further discussion below), then hab = gab + UaUb projects 
into the local rest spaces of comoving observers, where gat is the spacetime metric. The covariant 1 + 3 splitting of 
the Bianchi identities and the Ricci identity for incorporating Einstein's field equations as an algebraic definition 
of the Ricci tensor, Rab — Tab — ^Tgab, are the fundamental equations in the covariant perturbation approach. These 
equations may be written as propagation and constraint equations for covariant scalars, spatial vectors [Va — ha^Vb) 
and spatial 2-tensors which are symmetric and trace-free, i.e. which satisfy 

Sab = S(^ab) = ^a"^ hij"^ S (^^d) ^ \hcdS'^'^hab ■ 

Any spatial 2-tensor has the covariant irreducible decomposition 

Sab = \Shab + S(^ab) + ^abcS'^ , 

where S = h°-''Sab is the spatial trace and Sa = ^£abcS'"' is the spatial dual to the skew part. Here Sabc — Vabcdu'^ is the 
spatial permutation tensor defined by projecting the spacetime permutation tensor rjabcd- The covariant derivative Va 
splits into a covariant time derivative Aa-.- = u^S/b-A-a---, and a covariant spatial derivative Df,Aa... — h^'^ha'^ ■ ■ ■ Vd^c - - 



(Note that Dchab = = T>d£abc-) Then the covariant spatial divergence and curl are defined by [16 

div V = D''Va, curl Va = eafccD , 
(div^), = D'Sab , curl Sab = ecdia^'St)" . 

The fluid kinematics are described by the scalar = D"ua (expansion), the spatial vectors Ua (four-acceleration) 
and uJa — — ^curlua (vorticity), and the tensor <Jab = ^{a'^b} (shear). The locally free gravitational field is described 
by the electric and magnetic parts of the Weyl tensor, Eab — Cacbdu'^u'^ = E{ab) and Hab — \£ acdC'^'^ beu'^ = H/^ab)- 
The fiuid dynamics are given by the energy density p, the pressure p, and the dissipative quantities B (bulk viscous 
stress), Qa (heat flux, Qau"" = 0) and Hab = T^{ab) (shear viscous stress). These arise in the energy- momentum tensor 

Tab = PUaUb + (p + B)hab + '2q(aUb) + T^ab , (1) 

where 

qa = Qa+ (^) Ja (2) 



is the total energy flux relative to u"^, with n the particle number density and ja the particle flux {jau'^ = 0). The 
latter are combined in the particle four-flow vector 

AT" ^ nu'' + f . (3) 
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In a self-consistent thermo-hydrodynamic description, we need to introduce also the temperature T and specific 
entropy s per particle, defined in, or near to, equilibrium via the Gibbs equation 



T....(9+p.(i). (4) 



The hydrodynamic tensors Tab and N"" define two natural four-velocities - the particle (or Eckart) four-velocity 
Mp, for which — Q ^ Qa — qa, and the energy (or Landau-Lifshitz) four-velocity u^, for which ^ ^ Qa ~ 
— {p + p)ja/n. These four- velocities coincide in equilibrium, and differ by a small angle near to equilibrium.^ The 
four- velocity u°- may be chosen to be close to and Ug. Any small change in u"" produces second order changes 
(negligible in the linear regime) in p, p, n, T, and s |^. These scalars therefore coincide (to first order) with the 
corresponding scalars for a local equilibrium state. The bulk and shear viscous stresses B and nab are also invariant 
to first order under a small change in Both and j" undergo first-order changes, but the heat flux vector Q'^ is 
invariant to first order. From Eq. (0), we see that 

" " ^ q"- (5) 



° ip + p) 

For the isotropic and homogeneous Friedmann-Lcmaitrc- Robertson- Walker (FLRW) universes 

Da0 - BaP = DaP = D„B = D,n = BaT = D,s = , 

Ua = (^a ^ qa = ja = , 

Cfab = Eab = Hab = T^ab = . 

In covariant perturbation theory, a universe with small anisotropy and inhomogeneity is characterized by these quanti- 
ties being small, and one neglects terms which are nonlinear in them. Since these quantities vanish in the background, 
they are gauge-invariant |^. Note that FLRW models can admit scalar dissipation, in the form of a bulk viscous stress 
B (see, e.g., [p2|- pHT9[ | ) , reflecting the fact that expanding fluids in general cannot maintain equilibrium [Q. However, 
we shall follow the standard approach in irreversible thermodynamics of assuming an equilibrium background state, 
so that _B = in the background.^ For convenience, the linearized Bianchi and Ricci equations that underlie the 
covariant gauge-invariant theory are given in Appendix A (using the above notation and definitions, introduced in 
p6[ , which considerably simplify the original equations). Appendix A also contains useful differential identities. Note 
that in the background 

e ^SH, p^ 3H^{1 + K), H ^ ~\H^ [3(1 + w) + (1 + ?,w)K] ^ B^O^s, (6) 

where H — a/a is the Hubble rate, a is the cosmic scale factor, K = 0, ±(ai7)^^ is the dimensionless spatial curvature 
index, and w = p/p. 

Linearization of the number conservation equation Va-/V" = gives 

71 + 071= -D^ja . (7) 



Using the energy conservation equation (Al) and the number conservation equation (^, together with Eq. (||), the 
Gibbs equation (jj) implies that 

nTs = -3HB - D"Qa . (8) 

The contribution of shear viscous stress to entropy generation is via a nonlinear term (j°'^Tiab: so that in an almost- 
FLRW universe, the shear viscous stress does not contribute to s. Thus non-dissipative perturbations are not ade- 
quately characterized by s = 0. We need to specify that B = Qa = iTab = 0. 

Scalar perturbations are covariantly and gauge-invariantly characterized by the spatial gradients of scalars. Density 
inhomogeneities are described by the comoving fractional density gradient 



■^It has recently been argued [gyl that only the energy frame is suitable for the description of irreversible thermodynamics. 
^The consistency condition that this assumption imposes through the transport equation ( |4o| ) for B, is that the bulk viscosity 
should be much less than p9~^ . 
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= (9) 
P 

We define also the comoving expansion gradient and the dimensionless fractional number density gradient (not 
considered in |^,||), normalized pressure gradient, and normalized entropy gradient (see [Q) by 

Oa = aJJaO , = , Pa = , = • (10 

n p p 

Using the fact that p ~ p{p, s), and the Gibbs equation (^, we find 

Pa=C^Sa+rea, (11) 

ea ^ Sa ~ {1 + w)iya , (12) 



where the dimensionless quantities 



4-{^).r^±('^:). (13) 



dp J ' nT \ ds 



are respectively the adiabatic speed of sound and a non-barotropic index. Note that in equations (|T^) and (|T^), these 
quantities and w are evaluated in the background^ In the background 

c^ = ^, w^-3H{l + w){ci-w), (14) 
P 



where we have used s = and the energy conservation equation (Al). 

A covariant thermodynamic classification of scalar perturbations is as follows. Perturbations are non-dissipative 
if -B = Qa = TTab = 0, and then in particular s = 0, so that the specific entropy is constant along fluid flow-lines. 
If the specific entropy is the same universal constant along all flow-lines, i.e. if ea = in addition to s = 0, then 
the perturbations are isentropic, often (misleadingly) called 'adiabatic'. For isentropic perturbations, equations ( [Tl| ) 
and ( |l^ show that the number density perturbations and pressure perturbations are algebraically determined by the 
energy density perturbations: Va = '^a/(l +'w), Pa = c^Sa- The case of dissipative perturbations with Ca — 0, so 
that the speciflc entropy varies, but only along the fluid flow, will be called dissipative perturbations without entropy 
perturbations. The intcgrability condition Ca = 0, implies, via the gradient of the entropy evolution equation (H), that 



3HDaB + Da{B'Qb) ^0, (15) 

which is very restrictive, except in the case where only shear viscous stress is present. In general, dissipative pertur- 
bations will also involve entropy perturbations. 

The evolution of the temperature is clearly affected by the nature of the perturbations. In order to determine how 
this works, we use p and s as the independent thermodynamic variables in the Gibbs equation (^. The integr ability 
condition d^n/dpds = d^n/dsdp then gives 

where we used 

'dn\ _ _rfT_ 
.dsjp p + p' 

which follows from the Gibbs equation. Using the identity 

dT\ fdT\ 

s . 



dp ) ^ \9s J 



p 



*If B is nonzero in the background, then the background speed of sound acquires an additional dissipative contribution Cb, 
where (see [^) c^^ = /3o(p + p), and /3o arises in Eq. (^. 
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together with the energy conservation equation (Al) and the entropy evolution equation (^), the above equations lead 
to the temperature evolution equation 



T 
T 



1 



p{l + w) 



[3HB + D-^g,] 



nT^ V ds 



— ) [SHB + B'^Qa] 



(16) 



This equation reproduces the standard cooling rates for perfect fluids in the nonrelativistic and ultra-relativistic 
cases.^ In the general case, the source terms on the right hand side show the role of non-barotropic and dissipative 
effects. Note that the last term vanishes if the temperature is barotropic, i.e., if T = T(p). Bulk viscous perturbations 
counteract the cooling due to expansion, shear viscous perturbations do not affect the cooling rate (to first order), 
and the effect of heat flux depends on the sign of the divergence. For non-dissipative perturbations, the sign of the 
non-barotropic index r determines whether cooling is enhanced or retarded.^ 

We can also derive new evolution equations for the number density perturbations, entropy perturbations, and 
temperature perturbations. The com oving gradient of the number conservation equation (Q), together with the 
momentum conservation equation (A2) and the identity (A15), gives 



SrHi^a 



p(l 



3H DaB + qa+ 4Hqa + D°7r,f 



D,D^(Q6 



qb) 



(17) 



The comoving gradient of the entropy evolution equatio n (^ , together with the energy conservation equation (Al), the 
temperature evolution equation (16), and the identity ( A15 ), gives the evolution equation for entropy perturbations: 



e + 3H — w + r) e 



[3HD^B + (D^Qa) 



(18) 



where we have defined the scalar entropy perturbation 



e = aD^Ca = 



(19) 



Eq. dig) shows that for non-dissipative or shear viscous perturbations, entropy perturbations decay with expansion 
unlessc^ — w + r < 0. Defining the comoving fractional temperature gradient by 



oDgT 
T 



(20) 



we find from the evolution equation ( [l^ ) and the identity (A15), that the evolution equation of (covariant and gauge- 
invariant) temperature perturbations is given by 



-3 (cg + r) Haiia 
a id ^ - 



r) Oa - 3HaT)a (c^ 



pil 



w] 



3HDaB + Da D%fc 



r2 V alj^ 



3HDaB + T)a i^'Q 



(21) 



Now da contains more information than just the scalar density perturbations, since at each point, 5a picks out the 
direction of maximal inhomogeneity. The irreducible parts of the comoving gradient of Sa then describe completely 
and covariantly the variation in density inhomogeneities: 



aDtSa = hS hab + S.ab + SabcW 



(22) 



similar equation is given in but in the particle frame only, and without separating out the non-barotropic r terms. 
^Using the Gibbs equation (^, we can show that r = c'^[p{l -f w)a — ncp]/ncp, where a — n{dn~^ /dT)p > is the dilatation 
coefficient, and Cp = T{ds/dT)p > is the specific heat at constant pressure. 
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where the scalar part S = aD°'Sa = (aD)^p/p corresponds to the usual gauge-invariant density perturbation scalar 
£m 0, the vector part Wa = —^a curl (5a describes the rotational properties of inhomogeneous clustering, and the 
tensor part ^ab = aD{aSb) describes the volume-true distortion of inhomogeneous clustering. (These quantities were 
introduced in but only the scalar S was discussed.) These irreducible parts encode respectively the total scalar, 
vector and tensor contributions to density inhomogeneities. 

It is difficult to see how a rotation independent of the vorticity could arise, and indeed we can show that Wa is 
always proportional to the vorticity vector: 



Wa = ~3a^H{l + w)u;a ■ 



(23) 



This follows from the identity ( |A13| ) and the energy conservation equation (|Al| ). Thus rotation in clustering matter is 
inherited entirely (in the linear regime) from cosmic rotation: the vector part of density inhomogeneities is determined 
completely in direction by the cosmic vorticity. The expansion and pressure index affect the magnitude of the vector 
part. In particular, it follows that Wa — if the background is non-expanding or De Sitter (w = —1). 

The vorticity propagation equation ( A4) leads to the new evolution equation for the vector part of density inhomo- 
geneities: 



Wa + \H [3(1 - u;) -f (1 + ■iw)K] Wa 



curl 



qa + mqa 



(24) 



where we have used equations (g) and ( p^ , and the momentum conservation equation (A2) allowed us to evaluate 
curl Ma, together with the identity ( |A13| ). Unsurprisingly, Eq. (|^ shows that the scalar dissipative quantity B does 
not influence the evolution of the vector part of density inhomogeneities. |^ In the energy frame, heat flux also has no 
direct influence on Wa- Eq- ( p^ ) shows that for non- dissipative or only shear viscous perturbations, Wa decays with 
expansion unless 3(1 — w) + {1 + 3w)K < 0. For ordinary hydrodynamic matter, with < w < i, this inequality is 
never satisfied if the spatial curvature is non-negative. 

To derive evolution e qua tions for the tensor and scalar parts, we take the comoving spatial gradient of the energy 
conservation e quat io n (|AlD and the Raychaudhuri equation (A3), using the momentum conservation equation ( |A2| ), 
the identities ([A14|)- (|A16D , and Eq. (fl): 



6a^3wH6a-il+w)0a 



{1+W)9a 



-2H{l + w)ea - fi?^ [1 



qa + 4:Hqa 
>+ {l + w- 



p 



-DaDVb , 



(25) 



|c2) K] Sa - C^D^Sa 



-r [KH 



p 



[2,{l + w) + {l + ?.w)K] 



2p 



D TTafa 



Eq. ( psj ) can be shown to be in agreement with equation (61) of 
including bulk viscous effects. 

We can now decouple the equations: 

5a + H{2-Qw + icl) 5a - [l + 8w~ iw^ 



qa + ^Hqa 
cl curl Wa 
@, while Eq. (12 



(26) 

generalizes equation (62) of [b[ by 



clBHa 



-cl cuAWa + r{KH^ 



D2)e„ 



6c: 

a 



1 -3w 



'■ + lcl)K] Sa 



aH , r 
3 {qa^B^-iw 

p 



ici 



{l + iw)K] qa 

6H^ [1 - 3w -I- 2cl - (1 + 3w)K] qa - c^DaD^'gb} 

- |3i/D''7Tah + [2--3W + 3c2 - (1 + 3w)K] D^'TTat + DaD''DVb,} , 



(27) 



^Although the gradient of B occurs in the transport equation (^) for Qa, and therefore occurs on the right hand side of ( |24| ) 
in the particle frame {qa = Qa), the curl of this gradient is negligible by the identity (AIJ ), since B vanishes in the background. 
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where we have used equations (|), (|), (|ll|), and ([A1^)-(|A16|). The comoving gradient of the evolution equation ( |27| ) 
determines evolution equations for the scalar, vector and tensor parts of density inhomogeneities , incorporating all 
dissipative and entropy effects. We have alread y de ri ved th e vec tor e volution equation (^4|). Taking the comoving 
divergence of equation (^), and using identities ( Al4 )-( A16 ) and ( A20 ), gives the evolution equation for scalar density 
perturbations 



d + H{2-6w + 3cl) 6 - §i/2 [1 + 8w - ~ Gc^ + (l - iu? + 2(?^ K\ (5 - c^D^^ 
= S[e]+S[B]+S[9]+S[7r], 



(28) 



where the source terms arising respectively from entropy perturbations, bulk viscous stress, energy flux, and shear 
viscous stress are: 



S[e] = r (aXff^ + D^) e, 
S[S] = (3Xi/2 + D2)S, 

S[g] = 3ai? {g + H [1 - 9i(; + - (1 + 'iw)K\ q 

- f [1 + Bw - 9w^ - Bc^ + (1 - 9w^)K] q - c^D^} 
S[7r] = 3HS - [l + 6w - Sc^ + (1 + 3w)K] S + D^S , 

and we have defined the dimensionless perturbation scalars 



B 



S 



(29) 
(30) 

(31) 
(32) 



(33) 



Eq. (^8|) generalizes equation (74) of |^ to include bulk viscous effects, and is presented we believe in a more 
transparent form, which makes clear the physical meaning of each term. 

The new evolution equation for the tensor part of density inho moge ne ities follows from the trace-free symmetric 
part of the comoving gradient of Eq. (p7|), on using the identities ( A14| )-( A.16| ): 



U +{2-6w + 3c2) hU ~ IH^ [1 + Bw - 3^^ _ Gc^ + (l ~ 3w^ + fc^) K] U 

- c2D(<,Dt),5 = S[e]afc + S[B]ab + S[q]ab + S[7r]ah . 



The source terms are given by 

S[e]ab = SrKaH'^Di^aeb) 
S[BU 



2^2 

3K B,,Dh.B + DuB,.B, 



S[q]ab 



SM ab 



i 

P 



{Di^aiib) +H[l~2,w + icl - (1 + iw)K] D(,<7fc) 



[1 - 3u. + 2cl - (1 + 2,w)K] J^^aQb) ~ cl^ i^a^b)^" qc) , 

a2 

— {3iJD(aD=7Tf,)c + 3ij2 [2 - 3u; + 3c^ - (1 + 3w)K] D^^D'^tt,,), 



(34) 

(35) 
(36) 

(37) 
(38) 



Comparison of equations ( pS] ) and ( p^ shows that in the simplest case of isentropic perturbations, the density 
distortion tensor ^ab obeys the same equation as the scalar S, so that 



Cafc = Aab S , Aab = . 



(39) 



The presence of entropy or dissipative perturbations breaks the simple relation (|39|), and the evolution of the shape 
of density inhomogeneities is not directly determined by the scalar density perturbation. 



III. CAUSAL TRANSPORT EQUATIONS 



We are now ready to introduce the evolution equations obeyed by the dissipative quantities in the causal thermody- 
namics of Israel and Stewart M . This theory is based on a covariant treatment of the second law of thermodynamics 
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and the conservation equations, and its transport equations are confirmed by relativistic kinetic theory (via the 
relativistic generalization of the Grad approximation), which also provides explicit expressions for the various ther- 
modynamic parameters in the case of a dilute gas. The theory thus has a firm physical foundation. Furthermore, as 
pointed out earlier, dissipative signals propagate below the speed of light and the equilibrium states are stable, within 
the regime of validity of the theory. Thus the causal and stable thermodynamics of Israel and Stewart is a consistent 
relativistic thermodynamics which supercedes the non-causal and unstable theories first put forward by Eckart and 
Landau & Lifshitz. 

The predictions of the causal theory agree with those of the pathological theories in quasi-stationary situations. 
But when high-frequency/ short-wavelength effects are important, in the transient regime, the pathological theories 
are inapplicable. Thus these theories cannot cover the full range of behaviour of a relativistic fluid near equilibrium. 
Moreover, these theories cannot even constitute part of a consistent theoretical thermo- hydrodynamics because they 
are intrinsically not relativistic theories, given their pathologies. Thus our approach is to employ the causal thermo- 
dynamics to construct a self-consistent theory of cosmological density perturbations in the general case. In particular 
applications, where it can be argued that the non-causal theories will give reasonable results, we can then specialize 
the general equations appropriately. This is done in the companion paper H]. 

The full form of the transport equations, encompassing situations where the background equilibrium state is ac- 
celerating and rotating, and including terms which were neglected in the original theory and restored by Hiscock 
and Lindblom is given in Appendix B for convenience. Since we are dealing with cosmological perturbations, 
the background is non-rotating and non-accelerating, and spatial gradients of thermodynamic coefficients give rise to 
nonlinear terms. There are also linear terms, containing time derivatives of thermodynamic coefficients, which were 
restored by We will follow the arguments of which show that under many reasonable conditions, these 

terms may be neglected in comparison with t he o ther terms in the transport equations. n 



With these simplifications, equations (B1)-(B3) reduce to the causal transport equations 

(40) 
(41) 



B = -C[0 + l3oB - aoD'Qa 

Qa - ' 

TTab ^ -'2r][(Tab + P2'^ab ~ aiD(^aQb)] ■ (42) 



DaT + Tiia + Tf3iQa - TaoDaB - TaiB^iTab 



The coefficients C, k, and rj of bulk viscosity, thermal conductivity, and shear viscosity, appear also in the non-causal 
(and the non-relativistic) theories. The coefficients /3/ define characteristic relaxational time-scales 

To = CPo , n = K.Tf3l , T2 = 277^2 , 

which are often taken to be of the order of the mean collision time, but which are determined by coUisional integrals in 
kinetic theory Q. The non-causal theories are characterized by /?/ = 0. Intuitively, this corresponds to instantaneous 
relaxation to equilibrium when the dissipative 'force' is switched off. The coefficients a/, which also vanish in the 
non-causal case, arise from a coupling of viscous stress and heat flux (see Appendix B). They may also be found from 
kinetic theory in the case of a dilute gas. These transport equations hold in the energy and particle frames, with 
suitable simple changes in some of the thermodynamic coefficients (see Appendix B) 



The transport equations (40)-( |42|) are coupled to the evolution equations for density inhomogeneities - the scalar 
equation (p8|), the vector equation (^4D, and the tensor equation (^). They are also coupled to the evolution equations 
(p^, ( |l8| ) and ( pl| ) for number density, entropy, and temperature perturbations. In all of these evolution equations, 
except the entropy evolution equation (|l^), the energy flux vector qa occurs. Considerable simplification is thus 
achieved by choosing the energy frame {qa = 0), which is consistent with arguments in favour of that frame ^0 22 . 

In general, the coupling amongst the transport and evolution equations is highly complicated, although the coupled 
system can always be cast into a form suitable for numerical integration. Even in the non-dissipative case, when 
the transport equations fall away, the evolution equations themselves are coupled. In the simplest case of isentropic 
perturbations {B = Qa = iTab = Oj e = 0), the scalar equation ( p8| ) reduces to a wave equation only in 5 (with 
undamped phase speed Cg, as expected). In principle, the solution of this equation, and the solution Wa of Eq. ([24|), 
may be used to express Eq. ( |3^ ) as an equation only in ^ab- 

For non-dissipative entropy perturbations, decoupling S leads to a third-order equation. First we apply the operator 



3KH^ + to Eq. (Il8|), using identity (lAlTl) and 



Note liowever that it is not always reasonable to neglect these terms - see [[f2lll3] for examples 
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k = H(l + 3w) il + K)K, 



to get 



[{3KH^ + D^)e]' +H{2^3w + + 3r) [{3KH^ + D^) e] . 



Then we use Eq. ( p8| ) in this latter equation, together with identity ( |A17| ) and Eq. (j^) , to obtain the decoupled density 
perturbation evolution equation for non-dissipative entropy perturbations 



5 +[{A-%w + 6C3 +ir) H ^ {\nry]6 - \AH'5 - ^BH^S 



{3(c 



+ r)H~(\nry}c', + {ciy 



(43) 



where A and B are comphcated functions of w, c^, r, H, and K. In the case of a flat background K = 0, we have 
A = (1 + 84w - 27^2 - 69c2 + 27wc? - 12r + 36wr - 18rc? - 18cf ) H 



6(c^)' +2(2-6w-3c2) (Inr)' 
-1 + lOw - 39^2 + 24rw - ISc^ 4 
6 (c^)' - (1 + 8w - 3^2 ~ 6c^) (Inr)' 



B = (-1 + lOw - 39^2 + 24rw - ISc^ + 54wc,^ + 9w'^c^, + 3r - 9wV - 18rc^ 



18cf ) H 



A. The coupled system in general 

The complexity of Eq. (^3|) indicates the difficulty of decoupling the equations for dissipative perturbations. In 
general, the fact that the transport equations are first-order in time derivatives shows that any decoupling will produce 
at least one higher time derivative in the evolution equations. In the non-causal limit (/?/ = 0), when the derivative 
terms drop out of the transport equations, this does not hold, and the order of the equations is the same as in the 
non-dissipative case. 

The transport equations (^)-(^) contain further couplings amongst the dissipative quantities and couplings to 
the density and entropy and temperature perturbations. These further couplings are revealed when we take comoving 
spatial gradients, and use the following expressions: 



aD,C 



dp 



1 



1^ Sa + — 



nT \ds 



which follows from equation (|l 



3wHS -5 + 3HS 



which follows from Eq. (Eq), using the energy frame; 



1 



a(l + w) 



B + S] 



which follows from the momentum conservation equation (A2); 



which follows from the constraint equation ( |A6|) and the identity ( A20 ); and 



which follows from identity (AI8). Then operating on the transport equations (^)-(^2[) with, respectively, (a^/p)D^, 
(a/p)D", and {a^ / p)B°'B'' , we get the new causal transport equations for the scalar dissipative quantities 
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ToB + [1 - 3(1 + w)toH] B ^ {Caoa)I)^Q + 

Pii + w)[ \ 'c\dp 

tiQ+[1-3{1 + w)tiH]Q = 



3CH 



p{l + w) 
S 



S 

nT I as 



ap 



T 



kT 



ap{l 



[{aap{l + w)-l}B+ {aip{l + w) - 1} S] 



kT 



ap{l + w) 



\c^J + re] 



T2S 



4ry 

1 + w 



S - iwHS 



■i{l + w) 

We have defined the scalars for heat flux and temperature perturbations 



(44) 



(45) 



(46) 



Q = 



(47) 



The comoving spatial divergence of equation ( pT| ) gives the new evolution equation for scalar temperature perturbat 





f r\ r 


V 14 


- w ) . 



S + 3H{cl -w)S + 3Hre - SHa^B^ic^^ + r) 
„f^(S)[-«-°^e].3H(f^).. 



tons: 



(48) 



In summary, the coupled system that governs scalar dissipative perturbations in the general case is given by: the 
density perturbation equation (|2^), the entropy perturbation equation (|l8|), which we rewrite as 

(49) 



e + 3H(cl~w + r)e^ ~3HB - aD^Q , 



and the equations (^)-(Q) and (^). 

The number density perturbations do not occur in the coupled system. Once S, B, and Q are determined from 
the coupled system, the scalar number density perturbations, defined by v — aD^i^a, are found from the comoving 
divergence of equation (0). In the energy frame, this gives 



v + 'irHv= (l + w)~^ 



5 + ?,{ct~w + r)H5 + SHB + a^D^Q 



(50) 



The new evolution equation ( |50| ) shows how bulk viscous stress and heat flux govern the deviation of number density 
perturbations from energy density perturbations. Note that shear viscous stress does not directly afi^ect the number 
density perturbations. 

For specific applications, we present below the simplified coupled system that arises in two special cases when only 
one form of dissipation is present. 



B. Bulk viscous stress only 

The coupled system can be reduced to a pair of coupled equations in S (second-order in time) and e (second-order 
in time). In principle these may be decoupled. For a flat background, the equations are: 

d + H{2-6w + 3c,2) 6 ~ ^H^ (l + 8w - Sw^ - Gc^) 5 

= cjB^S+{w-cl)D'e^^D^e, (51) 

and 
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Toe + [1 - f (1 + 3w - 2c^ - 2r) tqH] e 



3H 



3H fdC 



r + 3(1 + w)rToH - tq [ci + r) + — ^ 

5. 



\ ^ 1 


i 3 


H(l + w) 


^+(l + u,) 



wQ + p{l + w) 



(52) 



Once these equations are solved for 5 and e, the other scalar quantities may be determined. Note that by the 
consistency condition (nS), the entropy perturbations cannot be zero unless B itself vanishes. 



C. Shear viscous stress only 



In the absence of entropy perturbations, the consistency condition (|l5| ) is identically satisfied if only shear viscous 
stress is present, and the system reduces to the pair of coupled equations: 



(2 - 6u> + 3<) 5 - [1 + Sw - 3w 



6c^ 



T2S 



l-2,w^ + 2cl) K] S 
c^D'^S + 3HS - [l + 6w- Sc^ + (1 + Sw)K] S + B^S , 

477 



1- H {'S{1 + w)t2- 



S 



4?7 



3(l + ii;) 



5 - 3wHS 



(53) 



(54) 



Finally, we point out an interesting feature of the case when only shear viscous stress is present (with or without 
entropy perturbations). The vector part Wa of density inhomogeneities satisfies the decoupled wave equation 



T2Wa + [1 + I {1 - W + (1 + 3w)K} T2H] Wa 



6(1 - w) - (1 + ly) (3 - w - &cl) T2H 



247? 



+ {2(1 + 3w) + [2(1 + &w + iw"^) - 6(1 + w)cl + 3(1 + iwfK] T2H 



81] 



1 - 



K 



Wa 



p{l + w) 



(55) 



The term D^ cTgh t hat arises from the transport e qua tion ( p^ ) is eliminated via the constraint equation (A6). We used 
the identity ( [A22| ) and the constraint equation (A7) to evaluate curl curl Wa- The undamped phase speed v is clearly 
given by 

,.2 V 



T2{P+P) 

It follows from the analysis of tej that w < 1, as expected from causality requirements. 
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APPENDIX A: COVARIANT PROPAGATION AND CONSTRAINT EQUATIONS 

The Ricci identity for and the Bianchi identities (incorporating the field equations via the Ricci tensor), may 
be covariantly spHt into propagation and constraint equations (see [D). In these equations are given in our 
streamUned notation in the exact nonUnear case, for a perfect fluid. For an almost-FLRW universe, with imperfect 
fluid, the hnearized form of the equations foUows. (In the non-dissipative case, all right hand sides are zero.) 



p+(p + p)0 = -3i/B-D%, (Al) 

{p + p)Ua + DaP = -D„B ~qa- AHqa ~ DV^fe , (A2) 

^+i6i2-D"w„ + i(p + 3p) = -fS, (A3) 

Ua + 2HLUa + jCUrlWa = , (A4) 

aab + '2Haab-'D{aUb)+Eab^ ^T^ab, (A5) 

lDad-D''(Jab+CmluJa=qa, (A6) 

D'^c^„ = 0, (A7) 

Hab - curlo-Qb - D(QCjf,) = , (A8) 

Eab + iHEab - CnvlHab + \{p+p)(Tab = -^T^ab - \H'Kab - {a%) , (A9) 

Hab + SHHab + curl Eab = IcurlTTah , (AlO) 

D^Eab - ^BaP = -Hqa - iO'^afc , (All) 

G'Hab - {p + p)uJa = -^cmlqa . (Al2) 

Some useful differential identities are P,p^ 

curlD,/ = -2/c^,, (A13) 

D2 (Daf) - D, (dV) + f (p - 3i/2) Baf + 2/curlw„ , (AM) 

(Da/)' = Baf-HBaf + fUa, (Al5) 

(aDaAb...)' ^aDaAb..., (A16) 

(dV)' = D^/ - 2i/DV + /D°?ia , (A17) 

D[,Dfc] V, = {H^ - ip) V[ahb]c , (A18) 

D[aDfc]5^'* = 2{H^- ip) Sia^^h]"^ , (A19) 

D°curlK = (A20) 

D^curl^ab = |curl (D^'Sab) , (A21) 

curl curl K = (oVb) - DV^ + 2 (ip - H^) Va , (A22) 

curl curl Sab = fD^.D^^,), - B^Sab + {p - 3H^) Sab , (A23) 



where the vectors and tensors vanish in the background, Sab — Si^ab) , and all the identities except ( |A13| ) are linearized. 
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APPENDIX B: FULL CAUSAL TRANSPORT EQUATIONS 



For completeness and convenience, we amalgamate the results of [Q [equations (7.1)]^ and [equations (21)- 
(23)][^ in our notation, and present the causal transport equations for viscous stress and heat flux in the general (near 
equilibrium) case, covering both cosmological and other (e.g., astrophysical) scenarios. 



B = -C\e + PoB - aoB^Qa 



(Bl) 



Qa = -kT —BaT + Ua+ PlQ (a) - a^^aB - aiD^TTat 



+ I (1 - 7o)™, (^) + (1 - 7i)r^/Dfc (^) - iQa 

TTab = -277 [o'ab + P2T^ {ab) - aiD(aQb)] 



(3ie + T 



T 



(B2) 



277 {a'lUl^aQb) + 2/32£cd(a7r6)'=w''} 



277 { liTQ(a'Db) ( Y j ~ 5'^"'' 



(B3) 



where = ha^Qh- The coefficients k, and rj are, respectively, the bulk viscosity, thermal conductivity, and shear 
viscosity. The relaxation coefficients /3o, and (32 are crucial to the causal behavior of the theory. The coefficients 
Qfo and ai arise from a coupling between viscous stress and heat ffux, as reffected in the entropy four-current [M 



1 

~ 2r 



1 

r 

1 



BQ"- 



T 

^PiQbQ' 



TTabQ 



ai 
Y 

- (il'^bc'^^'' 



Uqbu'^ 



2n-'qb) 



(B4) 



2T{p + p) 

The coefficients ao, Oq, ai, and a'l mediate the coupling of acceleration and vorticity to viscous stress and heat flux, 
while 7o and 71 appear due to a coupling of the spatial gradients of ai to viscous stress and heat flux. There are 
simple relations between the unprimed and primed a/ . 

A change from energy to particle frame results in a re-definition of various thermodynamic coefficients, but the 
transport equations maintain the same form. A partial comparison is given in [Q (p. 350), but the heat flux equation in 
the energy frame contains the spatial gradient of the thermal potential a = {p + p) /nT — s, and not the temperature 
gradient or acceleration, which arise in the particle frame. We can complete the comparison by using the Gibbs 
equation (m and the momentum conservation equation (A2) to show that 



nT 
P + P 



Baa 



1 

T' 



-B„T 



P + P 



B + B\r, 



(B5) 



in the energy frame. It follows that the energy-frame equation (2.38b) of [^ is in fact of the same form as the 
particle- frame equation (2.41b). 

In the cosmological setting, all the terms representing coupling of effects are nonlinear and may be neglected. To 
first order, Q(^a) — Qa and 'k(^ab) = '^ab- Furthermore, we follow the usual practice of neglecting the final term in 
square brackets on the right hand side of each transport equation. This can lead to problematic behavior in some 
cases, as shown in [|l2|,|l3| for the case of bulk viscosity. However, under a range of reasonable conditions, these terms 
may be neglected in co mpar ison with the remaining terms 1 14 2^ . With these assumptions, the terms in braces in 
the transport equations (B1)-(B2) all fall away, leading to the simplified cosmological transport equations (Eoh-(E2h. 



®Note that the ai in equation (7.1c) ofjfl] should be ai. 

^"Note that Va?" in equation (21) of pj should be q°'^Vaqb (using their notation). 
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